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As  shown  by  Pauli,  [Z.  Phys.  36,  336  (1926)],  the  electric  dipole  operator  r  can  be  replaced  by  the 
Runge-Lenz  vector  A  when  operating  within  the  n^  degenerate  manifold  of  hydrogenic  states  of  principal 
quantum  number  n.  We  seek  to  develop  similar  rules  for  higher  multipole  operators  by  expressing  equivalent 
operators  in  terms  only  of  the  two  vector  constants  of  motion — the  orbital  angular  momentum  L  and  the 
Runge-Lenz  vector  A — appropriate  to  the  degenerate  hydrogenic  shell.  Equivalence  of  two  operators  means 
here  that  they  yield  identical  matrix  elements  within  a  subspace  of  Hilbert  space  that  corresponds  to  fixed  n. 
Such  equivalent-operator  techniques  permit  direct  algebraic  calculation  of  perturbations  of  Rydberg  atoms  by 
external  fields  and  often  exact  analytical  results  for  transition  probabilities.  Explicit  expressions  for  equivalent 
quadrupole  and  octupole  operators  are  derived,  examples  are  provided,  and  general  aspects  of  the  problem  are 
discussed. 
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I.  INTRODUCTION 

Highly  excited  Rydberg  states  with  the  same  principal 
quantum  number  n  have  small  deviations  from  pure  hydro¬ 
genic  behavior.  The  degenerate  shell  of  these  states  forms  the 
basis  of  a  representation  of  the  0(4)  symmetry  group  [1] 
associated  with  the  Mr  Coulomb  potential  governing  the  dy¬ 
namics  of  the  Rydberg  electron.  Many  structural  properties 
of  the  Rydberg  atom  can  then  be  calculated  by  using  alge¬ 
braic  rules  and  group  representation  techniques.  These  fea¬ 
tures  combine  mathematical  beauty  with  pragmatic  useful¬ 
ness.  Moreover,  such  algebraic  techniques  facilitate  direct 
quantal  and  classical  solution  of  Rydberg  atoms  in  static  ex¬ 
ternal  electric  and  magnetic  fields  [2],  slow  collisions  with 
Rydberg  atoms  [3-5],  and  intrashell  dynamics  of  a  Rydberg 
atom  in  time-dependent  electric  and  magnetic  fields  [6].  For 
example,  analytical  probabilities  have  been  derived  [3-5], 
without  the  need  for  any  perturbative  and  numerical  analysis, 
for  the  full  array  of  transitions  in  atomic  hydrogen 

H{nl)  induced  by  a  time-varying  weak  electric  field  gener¬ 
ated  by  adiabatic  collision  with  slow  ions. 

The  dimension  of  the  degenerate  subspace  grows  as 
(without  electron  spin)  and  traditional  close-coupling 
(R-matrix)  calculations  using  spatial  wave  functions  become 
prohibitively  difficult  and  ultimately  impractical,  either  be¬ 
cause  of  the  sheer  dimension  of  the  space  or  because  of  the 
large  number  of  oscillations.  Rydberg  states  with  n  as  large 
as  several  hundred  are  now  accessible  to  observations  and 
experiments.  The  group  representation  technique  may  there¬ 
fore  offer  the  only  practical  and  effective  way  of  solving 
problems  involving  such  Rydberg  states.  In  so  doing,  some 
essential  underlying  physics  can  be  exposed,  as  an  additional 
asset. 

In  particular,  the  theory  of  a  Rydberg  atom  in  weak  exter¬ 
nal  electric  and  magnetic  fields  is  reduced  to  an  algebraic 


problem  that  is  extremely  well  poised  towards  extracting 
both  quantitative  analytical  results  and  qualitative  insight. 
This  is  possible  provided  the  necessary  operators  of  interac¬ 
tion  with  the  fields  can  be  represented  in  terms  of  integrals  of 
motion — the  angular  momentum  operator  L  (common  to  all 
central  field  problems)  and  the  Runge-Lenz  vector  A  (spe¬ 
cific  to  the  Coulomb  potential  alone).  This  key  step,  basic  to 
the  subsequent  algebraic  construction,  can  be  only  taken 
when  the  dynamics  becomes  restricted  to  the  subspace  of 
energy  degenerate  states,  i.e.,  when  the  external  perturbation 
is  so  weak  that  it  does  not  induce  n-n'  coupling  among  shells 
of  different  energy. 

All  of  the  above  work  [2-6]  was  based  on  a  paper  as  old 
as  quantum  mechanics  [7].  There,  Pauli  [7]  has  shown  that 
the  electric  dipole  operator  r  becomes  identical  with  the 
Runge-Lenz  vector  when  the  two  operators  are  restricted  to 
an  energy  shell  with  fixed  principal  quantum  number  n.  This 
can  most  easily  be  seen  by  comparing  all  the  matrix  elements 
of  these  operators  between  states  within  the  same  energy 
shell.  The  power  of  this  result  comes  from  its  general  valid¬ 
ity  and  utility  for  any  shell  with  quantum  number  n. 

The  advantage  of  expressing  the  intrashell  dynamics  of 
Rydberg  atoms  in  terms  of  the  A,L  set  of  constants  for  elec¬ 
tronic  Coulombic  motion  has  already  been  demonstrated 
[3-5]  for  collisional  /  mixing  transitions  induced  by  a  pro¬ 
jectile  charge-Rydberg  dipole  interaction.  Also  the  interac¬ 
tion  between  two  Stark-stretched  (polar)  Rydberg  atoms  has 
recently  been  expressed  [8]  in  terms  of  interactions  between 
the  permanent  multipoles  of  each  atom.  A  basic  question 
now  arises  quite  naturally  from  these  studies  [3-5,8].  It  is 
one  which  does  not  appear  to  have  been  previously  posed  or 
addressed.  Can  all  higher  multipole  interactions  be  equiva¬ 
lently  expressed  solely  in  terms  of  the  A ,  L  integrals  of  mo¬ 
tion  on  the  energy  shell? 
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In  an  effort  to  answer  this,  the  present  paper  considers 
general  equivalent  multipole  operators  in  Secs.  II  and  III.  A 
procedure  is  then  presented  and  applied,  with  examples,  to 
the  explicit  determination  in  Secs.  IV,  V,  and  VI  of  the 
equivalent  operators  for  the  dipole,  quadrupole,  and  octupole 
moments,  respectively.  The  algebraic  evaluation  of  the  op¬ 
erators  in  terms  of  (A,L)  is  conducted  in  the  Appendixes. 
Our  eventual  aim  is  to  provide,  if  possible,  the  full  solution 
for  general  multipoles.  Atomic  units  are  used  throughout  the 
paper,  unless  otherwise  noted. 


{ni,n2,in).  Stark  (n,k,m)  or  algebraic  {n,mi,m2)  quantum 
numbers  are  all  useful  hydrogenic  sets  (see,  e.g.,  Ref.  [8], 
Table  1).  When  an  operator  has  an  n-shell  equivalent  which 
can  be  expressed  in  terms  of  the  constants  of  motion  L  and 
A,  then  its  intrashell  matrix  elements  are  easily  calculated  in 
any  basis  of  states,  using  the  SO(4)  irreducible  matrix  repre¬ 
sentations.  For  example,  Pauli  [7]  has  proven  that 

3n 

=  - yA, 


II.  INTRASHELL  EQUIVALENT  OPERATORS 

The  angular  momentum  L  and  the  unrestricted  Runge- 
Lenz  vector 


U=^(pXL-LXp)-- 
2  r 

are  constants  of  motion  for  the  internal  Hamiltonian  H 
=p^l2-\/ r  of  the  Rydberg  atom.  A  more  convenient  form 
for  the  unrestricted  Runge-Lenz  vector  is 

U  =  -  p(r  •  p)  H-  rH. 


The  vector  operators  L  and  U  do  not  close  under  commuta¬ 
tion  to  form  a  Lie  algebra  because 


where  is  the  Levi-Civita  antisymmetric  symbol  for  indi¬ 
ces  ij  ,k=l  ,2,3.  If,  however,  the  action  of  operators  U  is 
restricted  to  the  Hilbert  subspace  of  states  with  principal 
quantum  number  n,  then  the  Runge-Lenz  vector  for  bound 
states  of  energy  £■„  can  be  defined  as 


A  = 


1  1  7  .  X 

==  -rp  -p{r-p)  +  rE„ 
V-  2£„  L  2 


(1) 


The  six  components  of  the  vector  operators  A  and  L  are 
generators  of  the  symmetry  group  SO(4)  of  proper  rotation 
in  four  dimensions.  They  satisfy  the  commutation  relations 
[Ai,Aj]  =  i€iji,L^,  [Li,Lj]  =  i€ijkLi,  and  [Li,Aj]  =  i€iji,Ak.  The  dis¬ 
crete  part  of  the  hydrogenic  spectrum  is  then  exhibited  [7]  by 
the  theory  of  irreducible  representations  of  SO(4). 

When  redefined  by  Eq.  (1),  the  Runge-Lenz  vector  A  acts 
only  on  states  within  the  n  shell  and  has  nonzero  matrix 
elements  only  between  states  within  the  n  shell.  If  P„  is  the 
projector  onto  the  n-shell  subspace  then  one  can  write 


A  =  P 


U 


"  J 


--P.,. 


V-  2H 


The  operators  A  and  U  /  V-2//  are  equivalent  because  all 
their  matrix  elements  are  equal  when  evaluated  between  all 
states  within  the  same  n  shell.  In  general,  the  operators  A 
and  B  are  equivalent  within  the  n  shell  if  all  their  intrashell 
matrix  elements  are  equal 

{ny\A\ny’)  =  {n'^B\ny') .  (2) 


The  quantum  numbers  y  label  the  basis  set  which  spans  the 
rP'  degenerate  subspace.  Spherical  {n,l,m),  parabolic 


so  that  the  dipole  operator  r  within  the  n  shell  is  equivalent 
to  -(3n/2)A.  The  dipole  matrix  elements  between  all  {l,m) 
states  of  the  same  n  are  therefore  simply  related  to  the  matrix 
elements  of  A,  which  are  then  algebraically  determined  most 
effectively  in  the  algebraic  or  Stark  bases. 


III.  MULTIPOLE  AND  MULTIPOLE-TYPE  OPERATORS 


Many  applications  require  calculation  of  matrix  elements 
of  higher  multipole  operators  and  therefore  it  is  useful  to  find 
their  n-shell  equivalents.  The  spherical-coordinate  represen¬ 
tation  of  the  multipole  operator  of  order  \  is 


Q 


477 

2\-h  1 


(3) 


where  r  is  the  electron  position  vector,  with  magnitude  r  and 
direction  f,  and  where  /jl  are  the  2X.H-1  components  with 
-X.<yU,<X.  Equation  (3)  for  X.=  l,2,3  provides  the  dipole, 
quadrupole,  and  octupole  operators,  respectively.  In  general, 
as  a  function  of  coordinates,  the  multipole  is  a  solution 
of  the  Laplace  equation  in  the  entire  free  space  (excluding 
the  singular  r=0  point). 

Another  definition  [9],  which  directly  reveals  the  irreduc¬ 
ible  tensor  properties  of  the  multipole  operators,  is 

•  -{{r  g>  r}®  ■  ■  ■  »  (4) 


which  represents  Q  as  a  multiple  irreducible  tensor  product 
of  vector  r  over  itself  X.  times.  A  multipole-type  operator  is 
obtained  if  the  identical  factors  r  are  replaced  by  different 
vectors  a,b,c,. . .  .  Theory  based  on  this  definition  will  not  be 
developed  here  but  will  be  the  subject  of  future  investigation. 

Representation  of  the  multipole  operator  by  its  Cartesian 
components  provides  several  advantages  for  the  approach 
taken  in  the  present  paper.  The  tensor  of  rank  X  with 
Cartesian  components 


(-  1 

^ - —d:  d:  A  - 

(2X-1)!!  ‘‘'2  '>^r 


(5) 


is  a  harmonic  polynomial  of  power  X  in  the  Cartesian  com¬ 
ponents  Xj  (j=l,2,3)  of  the  electronic  position  vector  r. 
Harmonic  polynomials,  by  definition,  satisfy  the  Laplace 
equation.  For  X  >  1 ,  the  tensor  (5)  can  be  expanded  as  a  sum 
of  terms,  as  a  result  of  taking  successive  derivatives.  One  of 
the  terms  is  x,-^, . . .  ,x,-^,  while  the  remaining  terms  contain 
at  least  one  &onecker  delta  symbol  for  a  pair  of  indices. 
However,  not  all  3^  components  of  the  tensor  are  inde- 
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pendent  because  the  tensor  is  fully  symmetrical  with  regard 
to  index  permutation  and  has  zero  trace  when  any  pair  of 
indices  is  contracted.  The  fully  symmetric  tensor  has 
independent  components  and  there  are  (j)  pairs  of  indices 
for  which  the  trace  is  zero.  The  tensor  has  therefore  only 
=  1  independent  components.  The  first  three 

multipole  operators  have  the  following  Cartesian  compo¬ 
nents: 


(2) 

4 


(3) 

Tljk  =  - 


Hi 

(6) 

1  2 

=  XiXj-y-8ij, 

(7) 

-r^{xiSj^  + XjSi,,  + , 

(8) 

for  the  dipole,  quadrupole,  and  octupole  moments,  respec¬ 
tively.  An  alternative  definition  for  the  Cartesian  components 
(5)  is  obtained  by  starting  from  the  monomial  ...  ,x,^ 

and  constructing  from  it  the  tensor  components  by 

adding  terms  such  that  the  result  has  both  the  required  sym¬ 
metry  and  the  zero  trace  condition  for  all  pairs  of  indices.  As 
an  example,  the  octupole  operator  is  explicitly  derived  via 
this  procedure  in  Sec.  VI.  This  is  also  the  way  one  can  con¬ 
struct  multipole-type  operators  starting  from  a  set  of  vectors 
(L  and  A  in  our  present  case)  which  replace  the  position 
vector.  For  example,  if  one  begins  with  vectors  a,b,c,..., 
the  multipole-type  operator  contains  the  Cartesian  compo¬ 
nents  bi^ ,Cj^, and  the  remaining  terms  are  obtained  by 
permutations  and  contractions.  Care  must,  however,  be  exer¬ 
cised  when  the  operators  do  not  commute. 

A  relation  between  the  spherical  (3)  and  Cartesian  (5) 
components  of  the  multipole  operator  is  facilitated  by  using 
the  following  definition  for  the  spherical  harmonics: 


2L-h  1 


4Tr(L  +  M)!(L-M)! 


1/2 


1 


SO  that 


r){L)  _ 

Cm  - 


(-  1)"(2L-  1)!! 


M 


M-k  k  L-M 


provides  the  required  result.  Explicit  relations  are 

for  the  dipole  operator. 


(9) 


Gi 


(2)  _  ±J2) 

0  -  ■ 


“(43 


G®  =  -  iqn  +  '2.iqn  "  >  (10) 

for  the  quadrupole  operator  and 

C?’  =  -  ^  +  Si,™  -  3,g!,  -  i,g,)  (11) 

for  the  octupole  operator.  We  now  seek  to  obtain  equivalent 
multipole  operators  within  the  n  shell  by  constructing  gen¬ 
eral  multipole-type  operators 

from  the  two  vector  constants  of  motion  L  and  A.  Because 
the  multipole  operator  of  order  X.  is  a  uniform  function 
of  coordinates  of  order  X,  coordinate  scaling  yr  implies 
that  the  operator  scales  as  We  therefore  re¬ 

quire  a  similar  property  for  the  equivalent  operator,  such  that 
T  contains  products  of  X  terms,  where  each  of  them  can  be 
either  L  or  A.  Further  restrictions  follow  from  parity 
(coordinate-inversion)  considerations — the  parity  of  is 
(-1)^,  the  parity  of  A  is  -1,  and  the  parity  of  L  is  +1. 
Multipoles  of  even  order  may  therefore  contain  only  prod¬ 
ucts  of  even  number  of  L  operators.  The  equivalent  dipole 
operator  (with  odd  parity)  is  expressed  only  in  terms  of  A. 
The  equivalent  quadrupole  operator  (even-parity)  has  terms 
A- A  and  L-L,  but  not  the  odd-parity  terms  with  A  L.  Be¬ 
cause  of  parity  considerations,  the  equivalent  octupole  opera¬ 
tor  can  only  contain  the  two  odd-parity  terms.  A- A- A  and 
A -L-L,  while  the  even-parity  term  A-A-L  is  forbidden. 


IV.  EQUIVALENT  DIPOLE  OPERATOR 

Using  the  general  procedure  outlined  in  the  previous  sec¬ 
tion,  we  seek  the  operator  equivalent  to  the  dipole  operator 
(6)  in  the  form 

=  PptPn  =  « A  +  (jL  . 

The  operator  L  has  the  opposite  parity  of  and  is  therefore 
precluded  by  setting  b  =  0.  The  coefficient  a  is  calculated  by 
comparing  the  matrix  elements  of  A  and  r.  Because  of  iden¬ 
tical  rotation  properties  of  the  two  vectors,  it  is  sufficient  to 
calculate  the  matrix  elements  only  along  one  direction.  It  is 
convenient  to  choose  this  direction  as  the  z  direction.  Be¬ 
cause  of  the  selection  rules,  the  dipole  operator  r  and  the 
Runge-Lenz  vector  A  have  nonzero  matrix  elements  only 
between  states  with  angular  momentum  quantum  numbers 
differing  by  one  unit  (/'  =  /±l). 
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Using  Eq.  (9),  the  intrashell  matrix  elements  of  the  z  com¬ 
ponent  of  the  dipole  are 


{nl  -  lm\z\nlm)  =  (f?„/-i|r|/?„/)(Y/_i„|Fio|F;J 

2”V  (2/-  1)(2/+  1)  ’ 


(12) 


where  the  results  (A3)  and  (All)  of  Appendix  A  have  been 
used. 

On  the  other  hand,  the  Runge-Lenz  component  has  the 
matrix  elements 


{nl  —  lm\Aj[nlm)  = 


jin^  -  f'){P'  -  m^) 

V  (2/-  1)(2/h-  1)  ’ 


(13) 


easily  deduced  from  Eq.  (A17). 

By  comparing  Eqs.  (12)  and  (13),  the  coefficient  is 
a  =  -3n/2.  This  reproduces,  as  expected,  the  result 
r^-(3n/2)A  originally  obtained  by  Pauli  [7].  The  position 
operator  r  can  be  therefore  replaced  by  the  Runge-Lenz  vec¬ 
tor,  when  restricted  to  states  within  the  n  shell.  This  result  is 
sometimes  referred  to  as  “Pauli’s  replacement  rule,”  which, 
in  addition  to  Pauli’s  original  proof  [7]  can  also  be  derived 
from  other  approaches,  as  in  Refs.  [2,3].  It  is  useful  to  note 
that  the  Cartesian  components  commute  when  unrestricted 
([jc,y]  =  0,  for  example),  but  behave  as  angular  momentum 
vectors  when  restricted  to  the  n  shell,  because  [A^,A^]  =  iLj.. 


V.  EQUIVALENT  QUADRUPOLE  OPERATOR 

The  quadrupole  operator  is  a  rank  2  tensor  with  Cartesian 
components  defined  by  Eq.  (7).  It  is  symmetrical  {q^- 
=  qi^i)  and  has  zero  trace  (S,q'„  =  0).  Two  quadrupolelike  op¬ 
erators,  symmetrical,  with  zero  trace  and  even  parity,  can  be 
constructed  from  A  and  L  and  can  contribute  to  the  equiva¬ 
lent  quadrupole  operator.  They  are 

0\]^  =  ^iAiAj  +  AjA,)-^A^Sij  (14) 

and 

0f  =  ^mLj  +  LjQ-^L^S^j,  (15) 

where  the  end  term  insures  zero  trace.  Mixed  terms  of  the 
form  AjLj  have  odd  parity  (sign  changes  under  coordinate 
inversion)  and  are  therefore  precluded.  The  equivalent  quad¬ 
rupole  operator  has  therefore  the  general  form 

P„q^^^P„  =  +  bO^^\ 

where  the  coefficients  a  and  b  are  determined  by  comparing 
the  matrix  elements  of  specific  tensor  components,  between 
states  with  angular  momentum  quantum  number  differing  by 
0  or  2  units,  i.e.,  (/'  =  /, Z±2). 

It  is  convenient  to  calculate  the  matrix  elements  of  the  zz, 
or  33,  components  of  the  quadrupole  and  quadrupolelike  op¬ 
erators  and  then  to  solve  the  set  of  equations 

{nlm\q^^^\nlm)  =  a(nZOT|(933^|n/m)  +  b{nlm\0^^j\nlm) , 


(n/'m|q'33^|nZm)  =  fl(n/'m|(?3*3^|n/OT)  +  b{nl'  m\0^^]\nlm) 
for  the  coefficients  a  and  b,  where  V  =  1-2. 


A.  Implementation 

(2) 

Explicit  calculations  of  the  above  matrix  elements  of  1733 , 
03*3^  and  (9 33^  yield  the  results  (A20),  (A21),  (A24),  and 
(A23),  derived  in  Appendix  A.  The  coefficients  are  thus  cal¬ 
culated  as  a  =  5n^/2  and  b  =  -n^/2.  The  quadrupole  equiva¬ 
lent  operator  can  therefore  be  expressed  exclusively  in  terms 
of  the  L  and  A  operators  as 


Pn^fPn  =  XiXj  -  -rSij  |P„ 


1.2 
'^J  3' 


2  ,  \ 

—  [AiAj+AjA,--A^S,j 


■J.  I 


+  (16) 


This  is  our  desired  rule  which  replaces  the  Cartesian  quad¬ 
rupole  components  (7)  by  Eq.  (16).  From  this  general  rule, 
the  equivalent-operator  relation 

=  y  (n^  +  3  +  Lj  H-  4A^  -  5A^) , 

(17) 

for  the  cylindrical  diagonal  element  squared  is  readily  de¬ 
duced.  Expression  (17),  originally  proven  by  Solov’ev  [10], 
was  key  to  theoretical  development  for  the  hydrogen  atom  in 
weak  magnetic  fields  and  in  crossed  electric  and  magnetic 
helds  (see  Ref.  [11],  for  example).  The  explicit  rule  (16) 
provides,  of  course,  all  five  independent  quadrupole  tensor 
elements. 


B.  Equivalent  operator  for  nlm—mlm'  transitions 


It  is  worth  noting  that  a  construction  was  developed  ear¬ 
lier  for  the  3D  rotational  SO(3)  group  with  generator  L  and 
the  unit  vector  n=r/r,  instead  of  the  SO(4)  group  vector  r. 
In  one  of  their  problems  [13],  Landau  and  Lifshitz  derive  the 
equivalent-operator  relation 


n.-ni.  —  Sjh  —  — 

'  *  3  (2/- 1)(2/h-3) 


Pi^k  +  Pk^i 


(18) 


to  be  compared  with  our  equivalent  operator  (16).  This  rela¬ 
tion  (18)  is  valid  for  the  Hilbert  space  of  states  on  a  unit 
sphere,  within  the  subspace  of  states  with  definite  total  an¬ 
gular  momentum  /  (but  with  different  m),  i.e.,  for  nlm 
^nlm'  transitions.  The  dehnition  of  equivalent  operators  for 
this  SO(3)  group  is 


{nl'}\A\nl'y'}  =  {nl'Y\B\nl'Y’}  (19) 

in  contrast  to  the  definition  (2)  of  hydrogenic  SO(4)  equiva¬ 
lent  operators.  The  present  study  provides  the  multipole  op¬ 
erators,  i.e.,  irreducible  tensors  built  from  electron  coordinate 
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r,  in  terms  of  the  hydrogenic  SO(4)  symmetry  group  genera¬ 
tors  L  and  A  which  permit  all  intrashell  nlm~^nl'm'  transi¬ 
tions. 

C.  Quadrupole  operator  in  the  SO(3)®SO(3)  representation 

In  the  Coulomb  problem,  it  is  conventional  to  replace  L 
and  A  by  a  pair  of  integrals  of  motion  and  related  by 

l(')  =  ^(L  +  A),  I®=^(L-A), 

l  =  i{i)  +  i(2)^  a  = 

The  operators  and  possess  all  the  properties  enjoyed 
by  independent  angular  momentum  quantum-mechanical  op¬ 
erators 


As  a  useful  application  of  the  present  equivalent-operator 
method,  we  evaluate  the  average  of  interaction  (20)  appro¬ 
priate  to  perturbation  theory.  Namely,  we  perform  in  the 
SO(3)  i8>SO(3)  representation  of  Sec.  V  C,  the  average  of  the 
quadrupole  operator  over  the  atomic  state  defined  by  two 
quantum  numbers  OTj  and  ni2,  where  m,,  1=1,2  are  the  re¬ 
spective  eigenvalues  for  projection  of  the  vector  operators 
onto  the  respective  axes  w,.  Let  a,-,  /  =  1 , 2  be  the  angle  be¬ 
tween  vectors  and  R,  and  let  (5  be  the  angle  between 
vectors  a>i  and  <«2-  Algebraic  manipulation  readily  yields  the 
important  result 


{mim2\qf]\mim2) 


"^2  2  2 
2mjcos  ofi  -  6m[m2COS  ajcos  a2  +  2m2COS  <2:2 


(I«)2  =  (i(2))2=j(,-+1),  j=^in-l), 

[l(l),l(2)]  =  0. 

These  pseudospin  operators  are  generators  of  SO(3) 
0  SO(3)  representation  of  algebra  of  SO(4)  symmetry  group 
in  quantum  mechanics  and  provide  the  algebraic  basis  set  [8] 
of  wave  functions  with  quantum  numbers  {n,mi,m2).  The 
equivalent  quadrupole  operator  rule  (16)  can  now  be  recast 
in  terms  of  pseudospin  vectors  and  as 

-  3n2[4»42)  +  /p)//)]  -  -  1)4 

+  2«2[i(1).i(2)]^.^^ 

D.  Application:  Averaged  ion-quadrupole  interaction 

In  a  number  of  physical  problems,  the  angular  momenta 
operators  and  are  quantized  with  respect  to  indepen¬ 
dent  axes  fe>i  and  W2  in  space.  One  example  is  the  hydrogen 
atom  in  crossed  electric  and  magnetic  helds,  where  the  vec¬ 
tors  ^a^  and  t02  are  expressed  [2]  in  terms  of  the  electric  field 
strength  and  the  magnetic  held  induction.  This  approach  is 
extended  also  to  time-dependent  helds  [6].  Another  example 
appears  in  the  theory  [3,5]  of  intrashell  mixing  in  excited 
hydrogen  atom  by  collision  with  a  particle  with  charge  Zg. 
Here,  in  the  corotating  frame,  the  electric  held  is  space-hxed 
and  directed  along  the  z  axis  towards  charge  Z^,  while  the 
effective  magnetic  held  is  perpendicular  to  it  and  normal  to 
the  collision  plane.  The  theory  [3,5]  of  intrashell  mixing  usu¬ 
ally  accounts  only  for  the  leading  and  dominant  charge- 
dipole  term  in  the  expansion  for  the  full  charge  (Z^) -Rydberg 
atom  interaction.  The  next  term  is  the  charge-quadrupole  in¬ 
teraction  which,  in  the  corotating  frame,  is 

3Zb^,  (20) 

where  R  is  a  vector  directed  from  the  atomic  nucleus  of  the 
target  atom  towards  the  projectile  of  charge  Zg  and  the  com¬ 
ponent  along  R  is  denoted  by  subscript  3. 


+  ((n^  -  l)/4  -  mjjsin^  +  {{n^  -  l)/4  -  /w^jsin^  q'2 


for  the  quadrupole  moment  averaged  over  an  algebraic  basis 
set  of  wave  functions  with  quantum  numbers  {n,mi,m2). 

VI.  EQUIVALENT  OCTUPOLE  OPERATOR 

The  Cartesian  components  of  the  octupole  moment  as  a 
tensor  of  rank  3  are  dehned  by  Eq.  (8).  The  octupolelike 
combinations  of  L  and  A  which  have  contributions  to  the 
equivalent  octupole  operator  within  the  n-shell  are  based  on 
A- A- A  and  L-L-A.  Other  combinations  have  even  parity 
and  are  forbidden  by  the  parity  rule,  as  previously  explained 
in  Sec.  III. 

Following  our  general  procedure,  we  seek  the  equivalent 
octupole  operators  in  the  form 

P„q^^'>P„  =  aO^'^'>  +  bO^^\ 

where  the  tensor  operator  is  derived  from  the  set 
A- A- A  and  O®  from  the  set  L-L-A,  respectively.  The  pre¬ 
cise  method  for  constructing  these  operators  is  now  ex¬ 
plained. 

The  matrix  elements  of  the  333  tensor  components  be¬ 
tween  states  with  one  and  three  units  difference  in  angular 
momentum  number  are  compared.  The  unknown  coefficients 
a  and  b  are  then  solutions  of  the  set  of  equations 

{nl—  lm\qf2i\nlm) 

=  a{nl  -  lm|6l3*33|nZOT)  +  b{nl  -  lm|C>33yn/m), 

{nl  —  3m\qf2^\nlm) 

=  a{nl  —  3m|(93*33|nZm)  +  b{nl  —  3m|C>333|«/m). 

The  Cartesian  components  O^jl  and  (9®  of  the  octupole¬ 
like  operators  are  now  constructed  from  AjAjA^.  and  LiLjA/^, 
respectively.  Three  operations  are  then  applied  to  these  el¬ 
ementary  combinations  to  insure  that  the  resulting  operators 
possess  the  following  properties:  (1)  they  have  zero  trace 
when  any  pair  of  indices  are  contracted,  (2)  they  are  fully 
symmetric  with  respect  to  index  permutations,  and  (3)  they 
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are  symmetric  over  the  order  of  noncommuting  factors. 

Given  three  vector  operators  a,  b,  and  c,  the  ijk  compo¬ 
nent  of  the  most  general  octupole-type  operator  is  given  by 


Oyi(a,b,c)  =  Z 


3  3 

1  1 
'Pfk~  <^ibsCs- 

J  .5=1  ^  i=l 


1 

5 


3 

^ik^J 

s=l 


(21) 


where  Z  represents  the  operation  of  symmetrization  with  re¬ 
spect  to  both  index  and  vector  permutations.  For  example, 

Z[aibjCk]  =  —2  2  pia)^{i)Pib).„(j)p{c)„(k),  (22) 

p  TT 

where  p  represents  a  permutation  of  the  vector  set  {a,b,c} 
and  where  tt  is  a  permutation  of  the  index  set  {i,j,k}.  It  is 
easy  to  check  that  the  combination  (21)  obtained  in  this  way 
possess  all  of  the  above  (l)-(3)  required  properties. 

Implementation.  According  to  the  general  prescription 
(21)  above,  the  operator  (8)  takes  the  form  Oy,(,=Oyi(r,r,r) 
so  that  its  zzz  (or  333)  component  in  Cartesian  representation 
is 


tlm  =  03nir,r,r)=.xl- 


3 

5 


as  expected.  From  vectors  L  and  A,  two  octupole  operators 
with  the  same  parity  as  for  can  be  constructed.  One  is 
0^jl  =  Oiji^iA,A,A)  and  the  other  is  (9®  =  Oy,(.(A,L,L). 

The  zzz  (or  333)  component  of  the  octupole-type  operator 
is 


^^3  ~  t?333(A,  A,  A)  —  Z 


^3  “  g  2  (2liA^)A3 


(23) 


With  the  aid  of  Eq.  (22)  and  the  commutation  relations 
[Ai,Aj}  =  iejji^Lp.  and  [Li,Aj}  =  i€jjpii^,  expression  (23)  reduces 
eventually  (see  Appendix  B  1)  to  simply 


tupole  operator  within  the  n  shell  is  therefore 


PnqfkPn 


35m 


15n^ 


Oijk{A,A,A)  +  Oy,t(L,L,A). 


(27) 


This  is  our  desired  rule  which  replaces  the  Cartesian  octu¬ 
pole  (8)  by  Eq.  (27),  where  the  general  components  are 
determined  in  the  prescribed  manner. 

In  particular,  the  sequences  involved  are  provided  in  Ap¬ 
pendix  B  for  the  0333  case,  as  an  example  characteristic  of 
the  overall  calculational  procedure.  Specifically,  Eq.  (27), 
with  the  aid  of  A^+L^+l=n^,  yields  the  expression 


Pnq'ilPn  = 


3L\A^  +  -jiA^A^+A^A^)  -  lAl 


■—{3n^  +  5) 


(28) 


for  the  equivalent  octupole  333-component  operator  in  the 
(L,A)  representation. 


VII.  LIST  OF  EQUIVALENT  MULTIPOLE  OPERATORS 
AND  A  TEST  EXAMPLE 

We  have  shown  that  the  Cartesian  dipole,  quadrupole,  and 
octupole  operators 

q^P  =  Xi,  (29) 

qf=XiXj-Y^Sij,  (30) 

qpl  =  XiXjXk  -  -r^Xjdjk  +  XjSii,  +  r^Sjj)  (31) 

have  the  equivalent  intrashell  dipole,  quadrupole,  and  octu¬ 
pole  operators 


0^33  =  ^3  -  ^  (^'^3  +  AaA^)  +  ^A3  .  (24) 

The  second  operator 
^333  ~  0333(L,L,  A) 

=  zi^LlA,  -  [(A^Qh  +  iLAs)L3  +  (44)A3]|  , 

(25) 

with  the  aid  of  similar  algebraic  reduction  and  the  additional 
identity  A-L=L-A=0,  eventually  reduces  (see  Appendix 
B  2)  to 

C®  =  lIa,  -  ^{L^A,  +  A,L^)  +  ^A3 .  (26) 

Using  Eqs.  (A26)-(A31)  of  Appendix  A,  the  solution  of 
the  set  of  equations  of  matrix  elements  yields  the  required 
coefficients  to  be  a  =  -^n  and  b=^n  .  The  equivalent  oc- 


qf^  =  -\nA„  (32) 

SA I  2  ,  \ 

+  (33) 

qpl  =  -J’^0‘^l  +  J'^0<§  (34) 

expressed  explicitly  in  terms  of  the  L  and  A  integrals  of  the 
motion  on  the  energy  shell.  Equation  (32)  reproduces  the 
Pauli  [7]  operator  replacement  rule  r  =  -(3«/2)A,  for  the 
dipole  Cartesian  component  (29).  Equations  (33)  and  (34) 
summarize  our  replacement  rules  for  both  the  quadrupole 
and  octupole  Cartesian  elements  (30)  and  (31).  The  operators 
0^^l  =  Ojjf.{A,A,A)  and  C>®  =  Oy,(.(L,L,A)  are  defined  in  Sec. 
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VI  and  are  calculated  in  Appendix  B.  The  scaling  law  is 
apparent  in  Eqs.  (32)-(34).  Note  that  equivalent  operators 
generally  contain  several  terms  with  simple  coefficients. 

Application'.  Multipoles  of  extreme  Stark  states.  As  a  test 
example,  consider  the  permanent  multipole  moments  associ¬ 
ated  with  the  extreme  Stark  hydrogenic  states,  i.e.,  for  those 
Stark  (parallel  and  antiparallel)  states  most  stretched  along 
the  positive  and  negative  direction  of  the  Z  axis.  With  the  aid 
of  the  replacement  rules  (32)-(34),  the  appropriate  spherical 
multipole  operators  are 

=  #  =  - yAs,  (35) 


r)(2)  -  ±M)  - . 
^0  “  2™  “ 


4  \  ^  3^ 


15«^ 


Ai  -  -A‘ 


(36) 


^(3)  _  fJ3)  _  _  . 


175m 


16 


Ai  +  -A,  -  — A^A,  -  — A,A' 


10 


10 


75«V  ,1  1  , 

+ -  L?A3  +  -A3  -  — L^At,  - 

16  \  ^  ^  5  ^  10  ^ 


—A.L^ 
10  ^ 


(37) 


The  parabolic  and  algebraic  representation  of  these  states 
are,  respectively, 

*/'(«-i)oo(r)  =  \  +  )  =  \n,j-j), 

<Ao(«-i)o(r)  =  \-)  =  \n,-jJ)- 

Any  of  these  hydrogenic  stretched  states  |  a)  =  |  ± )  has  expec¬ 
tation  value 

G!)'’  -  (38) 

For  the  “plus”  states  \+)  =  \nJ  ,-j)  of  the  algebraic  basis, 
we  have 


g^V)=(+ig;,‘'+>=-v^3(+), 


3(1) 


3n 


where  A3(+)  =  -2y=-(M- 1),  so  that 

G^'^(+)  = 


(1)  _  3n{n-l) 


is  the  permanent  dipole  moment  for  the  extreme  stretched 
Stark  state.  On  using  L3(+)=0,  L^(+)=m-1,  and  A\+)=n(n 
-1),  the  quadrupole  and  octupole  moments  are 

G®(+)  =  ^«V-1)(5«-7) 


and 


G^'’(+)  = 


l)(n-2)(7n-9). 


respectively.  All  these  moments  are  in  exact  agreement  with 
the  analytical  results  [8]  of  theory  recently  developed  solely 
for  the  case  of  these  extreme  Stark  states.  Matrix  elements  of 


operators  (35)-(37)  over  the  (n,  I,  m)  basis  have  also  been 
checked  for  n  =  1  -  6. 


VIII.  CONCLUSIONS 

We  have  shown  that  it  is  possible  to  construct  both  the 
quadrupole  and  octopole  operators  solely  in  terms  of  the 
A,L  operators  which  are  integrals  of  motion  on  the  energy 
shell  with  quantum  number  n.  We  have  provided  and  illus¬ 
trated  the  various  steps  involved  with  their  construction  and 
have  derived  explicit  expressions  for  these  operators.  The 
basic  importance  of  these  expressions  is  that  they  furnish  the 
ability  to  solve  various  structure  and  collision  problems 
solely  by  algebraic  group  theoretical  techniques  and  commu¬ 
tation  relations  based  on  the  SO(4)  symmetry  of  the  hydro¬ 
gen  atom,  when  the  dynamics  is  confined  to  the  energy  shell. 
We  gave  various  useful  averages  of  these  operators.  Our  gen¬ 
eral  theory  readily  provides  the  permanent  dipole,  quadru¬ 
pole,  and  octupole  of  polar  Rydberg  atoms  in  their  extreme 
Stark  states,  a  case  which  can  be  also  solved  by  less- 
sophisticated  standard  techniques.  The  present  treatment  will 
be  key  to  further  developments  in  the  theory  of  Stark  /  mix¬ 
ing  via  the  combined  charge-dipole,  charge-quadrupole,  and 
charge-octupole  interactions  evident  in  ion-Rydberg  atom 
collisions.  The  present  study  would  also  be  important  to  in¬ 
vestigation  [14,15]  of  Rydberg  atoms  in  the  field  of  electric 
multipoles. 

Although  the  dipole,  quadrupole,  and  octupole  are  the 
most  significant  in  many  practical  applications,  the  full  gen¬ 
eral  solution  for  any  multipole  remains,  at  present,  elusive. 
Although  plausible,  the  very  existence  of  equivalent  opera¬ 
tors  for  a  general  multipole  is  not  completely  certain.  Con¬ 
sider,  for  instance,  the  hexadecapole  operator  [i.e.,  Eq.  (5) 
with  X=4].  Here  three  tensor  operators  are  appropriate:  (5^^^ 
derived  from  the  set  A- A- A- A,  O®  derived  from  the  set 
A-A-L-L,  and  derived  from  the  set  L-L-L-L,  while 
the  sets  A-A-A-L  and  A-L-L-L  with  odd-parity  are  all 
forbidden.  The  equivalent  operator  can  then  be  taken  as  the 
linear  combination  aO^^^  +  +  cO^^\  Three  nonzero  non¬ 

diagonal  intrashell  matrix  elements  (/'  =  /,  V  =1-2,  and  /'  =  / 
-4)  are  then  used  to  fit  the  coefficients  a,b,c,  in  the  manner 
prescribed  in  Sec.  VI. 

Although  this  reasoning  could,  in  principle,  be  extended 
to  higher  multipoles,  the  overall  procedure  ultimately  be¬ 
comes  quite  cumbersome  for  actual  calculations  for  larger  X. 
The  simplicity  of  the  coefficients  in  formulas  (32)-(34)  pos¬ 
sibly  indicates  that  there  may  well  be  simpler  methods  of 
derivation.  Although  fruitful,  the  present  theoretical  system¬ 
atic  pole-by-pole  approach  may  not  be  sufficiently  powerful 
for  the  general  multipole  case.  Instead,  some  other  more  en¬ 
compassing  approach,  probably  based  on  study  of  commuta¬ 
tion  relations  (similar  in  spirit  to  the  derivation  of  Pauli  re¬ 
placement  rule  in  the  Appendix  of  Ref.  [3])  could  be 
developed. 

Finally,  we  indicate  that  there  is  another  context  where 
equivalent  operators  are  of  key  importance.  Namely,  higher 
order  contributions  from  external  fields  might  be  expressed 
in  terms  of  equivalent  operators.  These  effects  imply  virtual 
intershell  (^-changing)  transitions  conveniently  expressed 
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via  Green  functions.  Although  equivalent  operators  were 
constructed  quite  long  ago  by  Solov’ev  [12]  for  the  second- 
order  contributions  from  electric  fields,  higher  orders  have, 
as  yet,  not  been  considered.  Equivalent  operators  of  this  type 
are  beyond  the  scope  of  the  present  study. 
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APPENDIX  A:  INTRASHELL  MATRIX  ELEMENTS 


{nl\P\nl)  =  -n^{T>5n^  -  5n\6l(l  +  1)  -  5] 

8 

+  3/(/- l)(i +!)(/  + 2)},  (A7) 

{nl-  lir^|n/)=--n^(7«2-3/H5)\^^?^,  (A8) 

8 

{nl  -  2\r^\nl)  =  -n\ln^  -  /(/  —  1)  +  2] 

8 

XvK-/'V-(/-l)'],  (A9) 


{nl  -  3\r^\nl)  =  -  -  /2][„2  _  (/  _  \fjn^  _  (/  _  2)^] . 

(AlO) 


I.  Radial  and  angular  elements 

A  general  expression  for  intrashell  radial  matrix  elements 
of  integer  power  of  r  is  (see  Ref.  [16],  for  example) 


The  standard  angular  integral  is 


{nl'\r^\nt) 


^0 


Rnvir)Rni{r)r'^^dr 


■(«  +  />)!(«-/>- 1)!' 

U/ 

_(nH-/<)!(n-/<-  1)!_ 

1/2 


(8+1 


xEt 


(-  !)'■ 


X 


(«  +  /<  +  i)  !(«-/<-!  +  /) ! 


[n  +  />  +  i-(/3+  !)]![«-/>-  1 -Hi- (/3+  1)]!’ 


where  the  lower  limit  in  summation  is  io=max[0,/3+ 1 -(n 
-/>-!)],  />=max(/,/'),  and  /<=min(/,Z')-  This  general 
equation  yields  the  following  useful  matrix  elements  for  low 
powers  of  r. 


{nl\Y\nl)  =  1 , 

(Al) 

{nl\Y\nl)  =  -  /(/  -I-  1)], 

(A2) 

{nl  —  l|r^|/7/)  =  -  -nVn^  -  Y, 

(A3) 

(n/|r"|n/)  =  -n\5n^  -  3l{l  -I-  1)  -I-  1], 

(A4) 

{nl  -  l|r^|«/)  =  -  -n{5n^  -Y+  1)V«^  -  Y, 

(A5) 

{nl  -  2|r^|n/)  =  -n^V[n^  -  Y][n^  -  (/  -  1)^], 

(A6) 

/(2/i  +  l)(2/2+'Ty  ,0  ,,„3 

V  477(2/3-1-1)  /imj/jmj- 

From  standard  tables  (see  Ref.  [9],  for  example)  of  the 
Clebsch-Gordan  coefficients  ,  we  obtain  the  follow- 

1  j  /r*  j  1 2*' *2 ' 

ing  angular  integrals  for  the  quantum  numbers  of  interest: 

.  ,  ,  ^  I  ^ 

{I  -  lm\YiQ\lm)  =  \  — \- - — - -,  (All) 

'  ^  V  477  V  (2/-  1)(2/h-  1) 


(A13) 


/  I  I  ^  n~3 

(/  -  Im  T30  /m)  =  \l - - - — - - 

'  V  477  2  (2/- 3) (2/ H- 3) 


(2/-  1)(2/h-  1)’ 
(A14) 


(/  -  3m|T3o|/m) : 


4772(2/- 3)(2/-  1) 


X 


(/2  -  nYM  -  \Y  -  -  2Y  -  m^] 


{2l-5){2l+  1) 


(A15) 


2.  Dipole  matrix  elements  for  Sec.  IV 

On  using  the  relation  (9),  radial  (A3),  and  angular  (All) 
integrals,  the  dipole  matrix  elements  are 
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{nl  —  \m\q^^\nlm)  =  {nl  —  \m\Q^Q''\nlm) 


j  4-7T 

—  {nl  -  \\r\nl){l—  lm\Y i(\lm) 


3n 

2 


{n^-l^W-nY) 
(2/-I)(2/+I)  ■ 


(A16) 


3.  A  and  L  operators 


The  action  of  Runge-Lenz  vector  on  n-shell  states  is  (see 
Adams  [I],  p.  II2,  DeLange  and  Raab  [I],  p.  264,  or  Ref. 
[17],  for  example) 


A3|n/OT)  : 


(2/-  1)(2/+  1) 


\nl  -  Im) 


[n^-{l+iy][{l+lf-m^] 
{21+  l)(2/  +  3) 


\nl  +  Im). 

(A17) 


The  Pauli  replacement  rule  r-^-(3n/2)A  immediately  fol¬ 
lows  on  comparing  Eqs.  (A16)  and  (A17). 

In  terms  of  the  coefficient 


j{n^  -  yyP'  -  m^) 

V  (2/- 1)(2/h- 1)  ’ 


(A18) 


Eq.  (A17)  can  be  rewritten  as  the  linear  combination 

A-i\nlm)  =  g„i\nl  -  Im)  H-  g„,+i\nl  -I-  Im)  (A19) 

of  {n,l±l,m)  states.  The  operators  A+=Ai±L42  which 
change  m  by  ±1,  respectively,  can  also  be  written  as  the 
similar  combinations 


A+|n/m)  =  H-  I3,„  i_icjnl  -  Im  -I-  1)  -  r„,,;+iC„,+i|n/  H-  Im  H-  1), 

A_|n/m)  =  -  -  Im  -  1)  H-  r-„,,;+iC„/+i|n/  H-  Im  -  1) 

of  (M,Z±l,m±l)  states.  The  coefficients  are 

Pmj=  '{{l-m  +  l){l-m), 


7,71.1=  ^l{l  +  m+  l){l  +  m). 


c„,  =  v'(n2-/2)/(2/-  1)(2/h-  1). 

Eor  completeness  with  above,  the  components  L3  and  L+ 
=Li±iL2  of  the  L  operator  obey  the  standard  relations 

L3|M/m)  =  m\nlm). 


operators  and  frequent  use  is  made  of  the  standard 
commutation  relations  [A,-,Ay]  =  /eyj.Lj.,  and  \Li,Lj\  =  i€mLi^ 
together  with  [L;,Aj]  =  ieyj.A,(.  (which  relation  rotates  Aj  about 
axis  i  to  give  A]^,  thereby  confirming  the  vector  character  of 
A).  We  also  employ  the  additional  relations 
=  (A-L)  =  0  and  2^LjAj=  (L- A)  =  0  for  reduction. 

4.  Quadrupole  matrix  elements  for  Sec.  V 

On  using  the  relation  (10),  radial  (A4),  and  angular  (A12) 
integrals,  the  quadrupole  matrix  elements  are 

2 

(n/m|^®|n/m)  =  -(n/m|2®|n/m) 


2  A’Tt 

=  -  y  — (n/|r^|n/)(/m|F2olM 

/(/ H- 1)  -  3m2  2  .j,2  o,  13 

-(5n  -31  -31h-  1) 


3  {21-  l){2l  +  3) 

and  similarly,  using  Eqs.  (A6)  and  (A 1 3), 
{nl  —  2m|g33^|n/m) 


(A20) 


:  “  2m|Q®|n/m) 


2  477 


3  V  5 
2(2/-  1) 


{nl  —  2|r^|M/)(/  -  2m|  Y’2oI^^^*) 


X 


{n^  -  /2)[«2  _  (/  _  1)2](;2  _  _  ^2] 


(2/-3)(2/H- 1) 


(A21) 


On  using  Eq.  (A  19),  direct  algebraic  calculation  of  the 
matrix  elements  of  the  quadrupole  operators  (14)  and  (15) 
yields 


{nlm\0[l)nlm}  =  gl,  +  gl,^^  -  -  /(/  H-  1)  -  1], 


{nl  -  2m|e>^3V/'«)  =  gni8„i-i , 


(A22) 

(A23) 


(n/m|0®|n/m)  =  -  '^UO  +  1)  “  3m^],  (A24) 


L+lw/m)  =  w,„  /|n/m  +  1), 

L^\nlm)  =  u)^,„  i\nlm  -  1), 

where  &)„,/=  y(/-m)(/H-m-Hl). 

In  the  subsequent  reduction  of  the  basic  A/A,-  and  L,L, 

(1)  ^  (2) 

operations  within  the  quadrupole  operators  0]j  and  0]j  and 
the  AjAjAi^,  AiLjL),,  and  LjLjAi^  operations  within  the  octupole 


(«/ -  2m|(?33^|n/m)  =  0,  (A25) 

in  terms  of  the  g„i  coefficient  (A18).  The  relation  A^+L^ 
=  n^-l  has  been  used  for  Eq.  (A22). 

5.  Octupole  matrix  elements  for  Sec.  VI 

The  octupole  matrix  elements  calculated  using  Eqs.  (11), 
(A8),  (AlO),  (A14),  and  (A15)  are 
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{nl  -  Im|(7333|«/m)  =  “  \/  — {nl  —  l\r^\nl){l  -  Im|y3o|/m) 


3n^  (7«2  -  3/2  +  5)(/2  -  5m2  -  1)  (n^  -  P)iP  -  m^) 


8 


(2/-3)(2/  +  3) 


(2/-  1)(2/+  1) 


(A26) 


and 


{nl  —  3m|^333|M/m)  =  “  \/  — (nl  -  3|r^|«/)(/  —  3ot|  Y30I//W) 


/4t7 


35n^ 


8(2/-3)(2/-  I) 


(^2  -  /2)[„2  _  (/  _  1)2][„2  _  _  2)2](/2  _  _  2)2  _  ^2] 

(2/-5)(2/+I) 


(A27) 


Direct  calculation  based  on  Eq.  (A  19)  then  provides  the 
following  octupole  matrix  elements: 


{nl  —  lm|(93*33|n/m) 

“  8nl 


sh-i + sh + sIm  +  ^  ^  - 1) 


(A28) 


{nl  -  3m\0%\nlm}  =  gnigni-\gni-2^  (A29) 


{nl  —  lm|(9333|n/m)  =  -  ~^gnM^  ~  5^^  “  1)>  (A30) 

(«/ -  3m|(93|3|n/wr)  =  0,  (A31) 

in  terms  of  the  coefficient  (A18).  The  relation  A^+L^=n^ 
- 1  has  been  used  in  Eq.  (A28). 


APPENDIX  B:  (A,L)  Representation  for  operators 

^333  ^333 

We  illustrate  here  the  procedure  for  evaluating,  in  terms 
of  ( A ,  L) ,  the  03*33  and  6)333  operators  from  their  basic  defi¬ 
nitions  (23)  and  (25).  At  the  outset,  proper  account  must  be 
taken  of  the  inherent  2^-operation  prescribed  by  formula  (22) 
for  symmetrization  with  respect  to  both  index  and  vector 
permutations. 


Z'^  (A^JA3  -  2  (AiA^3  +  AjA3Aj  +  A^/l3  H- A  jA3Aj 

.  6  , 

+  A3A5A^  +  A3A^j) 

=  ^(A^Aj  -I- A3A2)  H-  A/I3A, 

=  -(A2A3  +  A3A2)  +  “2  (AjA^3  -  A^5A3 

3  6  , 

+  2A^3Aj  -  +  A3AjAj) 

=  -(A2A3  -I- A3A2)  +  “2  (Aj[A3,AJ 
2  6  , 

-  [A3,A  JA^) 

=  -(A^A^  -I- A3A2)  +  -[A1/L2  -  /L2A1 
2  6 

+  A2(-  ILi)  -  (-  /Li)A2] 

=  -(A2A3  -I- A3A2)  +  -/([Ai,L2]  -  [A2,Ti]) 

2  6 

=  -(A2A3  -H  A3A2)  +  -i{iA^  +  lA^) 

2  6 

=  -(A2A3-I-A3A2)--A3. 

2  3 

The  6)3*33  operator  (Bl)  therefore  reduces  to 


1.  The  6)333  component 
The  definition  of  6)3*33  is 


6)3*33  -  A3  -  ^^(A^A3  +  A3A^)  +  ^A3 


(B2) 


C'^33  =  0333(A,A,A)  =  Z 


3 

“2  (AjAjA3 

c 


(Bl) 


The  hrst  term  ZA3=A3  is  left  unaffected.  The  symmetrized 
sum  can  be  evaluated  via  the  following  progression  of  steps: 


as  stated  by  Eq.  (24)  of  the  text. 

(2) 

2.  The  6)333  component 

(2) 

6)333  is  defined  by 
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C>333  -  0333(L,L,A)  —  zi  [.{AsLs)L^  +  (Ly4j)L3  +  (LjLj)A3]  i . 


(B3) 


Because  the  operators  L3  and  A3  commute,  the  hrst  term  ZL\A^=L\A^  is  left  unaffected.  Symmetrization  of  the  remaining 
terms  proceeds  via  the  following  steps: 

-22  +  (■^-i^s)^3]  -7^  i^s^sL^  +  A^L^L^  +  ijA^Z,3  +  L^L^A^  +  L^A^L^  +  +  L^L^A^  +  A^L^L^ 

s  6  , 

+  +  L^A^Lg  +  L^L^A^  +  ^3^s^s) 

=  “2  (^jf-if-3  +  ^i-^-3-^-i  +  ^As^3  +  L.AAs  +  +  ^3^i^i)  +  3  +  ^3^^) 


+  2 2  ^A3^s ~  2^  +  L..L2A3  +  ^ {iAa^  +  A^A)  +  2 ^  ■ 


(B4) 


The  hrst  sum  on  the  right-hand  side  of  (B4)  transforms 
from 

51  =  2  AsA^s  +  As^3^s) 

s 

=  2  AsA^s  ~  ^s^sA  + 

s 

=  2(A,[L3,4]-[L3,LJA,) 


■  +AA)  +  ^{LiiA2  —  iA2A  ~  Ai^\  + 

■  ~iAA2  +A2A)  +  “(i[Ti,A2]  -  i\L2A\\} 

:  ^(l2A3  +  A3L2)  +  i(/iA3  +  (-  0(-  iW  (B5) 


to 

5]  =  A]iL2  -  iL2Ai  H- A2(—  Of-i  “  (~  Of-1^2 

=  i[Ai,L2]  ~  *[-^25^1]  =  “  2A3. 

The  second  sum  on  the  right-hand  side  of  (B4)  becomes 
reduced  by  the  steps  from 

*^2  =  2  l^a^A 

s 

=  “2  AsLA3  ~  i-AA3  +  2LjA3Lj  H- A3L^Lj  -  A3LjLj) 


^{AA^  +A2A)  +  2^  (-^-i[-43>f-j]  “  [-AsjTJLj) 


to 


S2=-(L^A3+A3L2)-A3. 


The  full  O333  operator  (B3)  is  therefore 


1 


0%  =  LiA2-Y^{LA2+A2A)  +  -A 


(B6) 


as  stated  by  Eq.  (26)  of  the  text.  The  end  term  of  Eq.  (B6) 
being  extracted  from  the  detailed  algebraic  calculation 
above,  is  the  least  obvious. 
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